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The O(N) model of layered antiferro- and ferromagnets 
with a weak interlayer coupling and/or easy-axis anisotropy 
is considered. A renormalization group (RG) analysis in this 
model is performed, the results for N = 3 being expected to 
agree with those of the 1/M expansion in the CP M ~ 1 model 
at M — 2. The quantum and classical cases are considered. A 
crossover from an isotropic 2D-like to 3D Heisenberg (or 2D 
Ising) regime is investigated within the 1/N expansion. An- 
alytical results for the temperature dependence of the (sub- 
lattice) magnetization are obtained in different regimes. The 
RG results for the ordering temperature are derived. In the 
quantum case they coincide with the corresponding results of 
the 1/N expansion. The numerical calculations on the base 
of the equations obtained yield a good agreement with ex- 
perimental data on the layered perovskites La2Cu04, K2NIF4 
and Rb2NiF4, and the Monte-Carlo results for the anisotropic 
classical systems. 



I. INTRODUCTION 

The problem of layered magnetic systems is of inter- 
est both from theoretical and practical point of view. 
Here belong, e.g., quasi-two-dimensional (quasi-2D) per- 
ovskites , ferromagnetic monolayers and ultrathin films 
|^| . Such systems possess magnetic transition tempera- 
tures which are low in comparison with the intralayer 
exchange parameter J and are determined by magnetic 
anisotropy and/or interlayer coupling. 

The crucial role in the systems with small interlayer 
coupling (or anisotropy) belongs to the temperature 
crossover from a "2D-like" (isotropic) regime to the crit- 
ical 3D (or 2D Ising) regime respectively (see, e.g., Refs. 

In the general case where both interlayer cou- 
pling and anisotropy parameter are present, the situa- 
tion is more complicated: with increasing temperature 
the 2D-like Heisenberg behavior changes to the 2D-Ising 
or 3D-Heisenberg one depending on that anisotropy or 
interlayer coupling dominate, and finally system passes 
to the 3D Ising behavior. 

There exist a number of approximations which treat 
thermodynamics of layered systems. The standard spin- 
wave theory (SWT) describes satisfactorily the region of 
rather low temperatures only. Somewhat better results 
can be obtained by taking into account the temperature 
renormalization of the interlayer coupling parameter and 
the anisotropy parameter. The temperature dependence 
of the anisotropy parameter within the spin-wave theory 



was considered in Refs. [^-j6[. More systematic way to 
consider such renormalizations is the self-consistent spin- 
wave theory (SSWT) 0,§ which was applied to quasi-2D 
and anisotropic 2D magnets in Refs. 

SSWT takes into account the interaction between spin 
waves in the lowest Born approximation. However, at not 
too low temperatures this approximation is insufficient. 
In particular, the values of the ordering temperature in 
SSWT are still too high in comparison with experimental 
ones, and the critical behavior is quite incorrect. Thus 
the summation of leading contributions in all orders of 
perturbation theory should be performed. At the same 
time, to describe the behavior of the order parameter in 
the critical region one has to take into consideration fluc- 
tuation (non-spin-wave) contribution to thermodynamic 
quantities. Ising-like excitations in a classical anisotropic 
model were considered in Ref. |u|. At the same time, in 
the quantum anisotropic case this treatment meets with 
difficulties |TJJ|, and 3D fluctuations in the critical re- 
gion of quasi-2D magnets cannot be considered in the 
approach j[3 14 too. 

A possible way to sum up an infinite sequence of 
perturbation contributions is the renormalization group 
(RG) analysis. The RG approach was successfully used to 
consider the classical isotropic magnets with the space di- 
mensionality d = 2+e [ plpl . In this case the fixed point 
is small, T* ~ eJ (J is the exchange parameter), and the 
standard technique of the e-expansion can be applied. 
Physically, this means that the excitation-spectrum pic- 
ture differs somewhat from the spin-wave one (as dis- 
cussed in Ref. ||, the fluctuation corrections to the ex- 
citation spectrum read SEq ~ e In q) . The RG method 
was applied also to quantum 2D isotropic magnets in Ref. 

The scaling behavior in the quasi-2D or anisotropic 2D 
systems is expected to differ from the isotropic d = 2 + e 
magnet. In these cases the renormalized value of T* / J at 
the stability point of the RG transformation is not small, 
which corresponds to above-discussed crossover from the 
2D-like Heisenberg to 3D Heisenberg (or 2D Ising) critical 
regime. Thus the RG method does not work when pass- 
ing to the true critical region. The latter region should to 
be considered with account of essentially non-spin-wave 
(fluctuation) excitations. 

To take into account non-spin-wave (fluctuation) ef- 
fects it is convenient to use, instead of the original Heisen- 
berg model, models with large degeneracy, which enables 
one to introduce a formall small parameter in the the- 
ory. In Ref. M an isotropic quantum quasi-2D antifer- 
romagnet was investigated within the 1/N expansion in 
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the quantum O(N) model pS] (in the Heisenberg model, 
N = 3). It was demonstrated that the renormalizations 
of the interlayer coupling parameter in comparison with 
the usual spin- wave theory (or SSWT) determine consid- 
erable lowering of the transition point. The same situa- 
tion should be expected for 2D anisotropic magnets. At 
the same time, the 1/N expansion meets with some dif- 
ficulties at the description of the 2D-like fluctuations in 
the renormalized-classical temperature region near d = 2 
where the series expansion is performed in powers of 
1/(N - 2) rather than 1/N @|. 

Thus the RG approach and 1/N expansion in the O(N) 
model are expected to have advantages in different tem- 
perature regions. Whereas the first method describes well 
the 2D-like regime, the 1/N expansion treats satisfacto- 
rily the critical region. An advantage of the RG approach 
in comparison with the technique of the 1 /N expansion 
in the O(N) model is that it enables one to consider the 
quantum ferromagnet case where the partition function 
cannot be generalized to arbitrary N. The RG analysis 
permits also to treat the classical case, which is difficult 
within the 1/N expansion. (In the classical case, there is 
no natural upper cutoff for quasimomenta, which is the 
temperature in the quantum case, and the original lattice 
version of the partition function should be considered.) 

The Heisenberg model can be also considered as a par- 
ticular case of the SU(M) model (or of its continual ana- 
log, CP M - X model) with M = 2 (see, e.g., @). Since 
the M — > oo limit corresponds to SSWT (see, e.g., Ref. 
H), at finite M thermodynamics is described in terms 
of spin-wave picture of the excitation spectrum. The 
corresponding 1/M-expansion contains in the 2D case 
infrared-divergent terms [ pCfl and is rather inconvenient. 
For d not close to 2, this expansion does not enable one to 
obtain correct values of the critical exponents [^J. Thus, 
as well as in RG approach, only the 2D-like isotropic 
Heisenberg region can be considered within the 1/M- 
expansion. Furthermore, the results of CP M_1 model 
to n-th order in 1/M are expected to coincide at M = 2 
with the (n + l)-loop RG analysis for 0(3) or, equiva- 
lently, CP 1 model (we do not know a general proof of 
this statement, but this is true in the d — 2 + e case for 
n = 0,1, see Ref. fl)). 

The simplest one-loop RG analysis was applied earlier 
to calculation of the Curie temperature of anisotropic 2D 
ferromagnets p2fl . As we shall see below, the results of 
this work become considerably modified by the two-loop 
corrections which were not taken into account in Ref. 
p2j. The momentum-shell version of the one-loop RG 
approach in the classical quasi-2D magnets was consid- 
ered in . However, the authors of this paper passed to 
continual limit in the direction perpendicular to layers, 
so that results of Ref. [^3| at small interlayer coupling 
have qualitative rather quantitative character. 

In the present paper we consider thermodynamics of 
the quantum and classical layered magnets with small 
interlayer coupling and anisotropy within the consistent 
two- loop RG approach and first-order 1/N expansion in 



the O(N) model. 

The plan of the paper is as follows. In Sect. 2 we con- 
sider the continual and lattice versions of the 0(N) model 
for anisotropic layered ferro- and antiferromagnets. In 
Sect. 3 we apply the renormalization group approach 
|l5| , ^7| to the quantum O(N) model in the (isotropic) 2D- 
like region (the classical case is considered in Appendix 
B). In Sect. 4 we investigate the same problem within the 
1/N expansion in the O(N) model up to the first-order 
and treat the crossover from the 2D-like to 3D Heisen- 
berg (or 2D Ising) behavior. In Sect. 5 we summarize our 
results and compare them with experimental data on lay- 
ered antiferromagnets. 



II. CONTINUAL AND LATTICE MODELS FOR 
THE SPIN SYSTEM 

We consider the Heisenberg model with small inter- 
layer coupling and easy-axis anisotropy 
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where J > for a ferromagnet, J < for an antifer- 
romagnet, <$ii and 5± denote nearest neighbors within a 
layer and for different layers, a > is the interlayer cou- 
pling parameter, 77, £ > are the parameters of the ex- 
change anisotropy and single-site anisotropy respectively. 
The partition function of the model (Q) can be repre- 
sented in terms of a path integral over coherent states 
(see, e.g., Refs. @fl): 



DnDAexp 



JS 2 



l/T 



2 i , , d\\i . 



+ -n t n l+s± + r?nfn| +5|| + sgn(J)C«)^ + hn\ + i X^nf - 1) 

with nj(r) the three-component unit-length vector field, 
A(n) the vector potential of the unit magnetic monopolc, 
which satisfies the equation V x A(n) • n = 1, £ = 2£(1 — 
1/2S) and the summation over <5ii,<5j_ in (Q) is assumed. 
We have also introduced in (Q) the external magnetic field 
h to perform the calculation of spin correlation functions. 
The term with the time derivative corresponds to the 
Berry-phase contribution pq| . Depending on the value 
of T/JS, two cases are possible: (a) the classical case 
T > JS and (b) the quantum case T < JS. 
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Consider first the classical case. The main contribu- 
tion to (^J) comes from time-independent paths, and the 
partition function reduces to 



DnDX exp 




(5) 



where <Tj a is the three-component unit-length field and 
c = \/8JS is the bare spin-wave velocity. Unlike the 
quantum ferromagnet case, this model can be trivially 
extended to the O(N) symmetry with arbitrary TV, the 
iV-component vector field <7j = {<ji...(tn} being intro- 
duced and a z being replaced by <tn- Writing down 



+VnM +h + C«) 2 + hn\ + i A(n 2 - f )] } 

with p° = | J | 5 2 the bare spin stiffness. To derive 
in the antiferromagnetic case we have replaced 



a = {m. 
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-Xi at one of two sublattices. Thus in the 



classical case the results for Z are identical for a ferro- 
and antifcrromagnet. In the continual limit the parti- 
tion function (|J) coincides with that of the well-known 
classical nonlinear-sigma model ||. However, if one is 
interested in thermodynamics in a wide temperature in- 
terval (not only in the critical region) , the continual limit 
cannot be used since not only long-wave excitations con- 
tribute thermodynamic properties. 

Now we treat the quantum case. Then the temper- 
ature plays the role of a natural upper limit cutoff for 
frequencies of the fluctuations. Thus we may pass to the 
continual limit within each layer. For a ferromagnet we 
use the representation 
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where <9 A 
shift iX - 



= (d/d(cor), V), and we have performed the 
>iX + f. 



A(n) 



1 + (zn) 

(z is the unit vector along the z-axis). Then we obtain 

l/T 

i /<:.: 

: exp 
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III. RENORMALIZATION GROUP ANALYSIS IN 
THE 2D-LIKE REGIME FOR THE QUANTUM 
CASE 

Using the above expressions for the partition function 
we can develop a scaling approach. We use the field the- 
ory formulation of RG transformation |27||l5| ] . To develop 
this approach we pass to the renormalized theory with 
the use of the relations |l5|j 
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.9 = 9rZ\, u = ub.Z u 
n = n R Z, h = h R Z 1 / 
f = IrZ2, a = a R Z 3 



Z 
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where the index R corresponds to the renormalized quan- 
j^stities, the bare coupling constant g and the dimensionless 



where i z is the number of a layer, 7r 
two-component vector field, 



n — (nz)z is the 



inverse temperature u are determined by 

1,1 g = 1/5, u = JS/T (FM) 
g = c /p°, u = c /T (AFM), 



(13) 
(14) 



/ = C + 477 = 2C(1 - 1/25) + 4r? 



(8) 



and we have made the shift iX — ► iX+f before integrating 
over A. 

In the antiferromagnetic quantum case we use the Hal- 
dane mapping J2f| (see also Ref. || ) to integrate over the 
"fast" components of n. Thus we pass to the partition 
function of the quantum nonlinear sigma model 
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The renormalization constants Zi are chosen from the 
condition that the thermodynamic quantities are inde- 
pendent of an upper cutoff. Since the nonlinear-sigma 
model is renormalizable (see Ref. |27j]), five renormal- 
ization constants for five independent parameters of the 
model are sufficient to this end. To calculate renormal- 
ization constants it is sufficient to calculate the renormal- 
ization of the o ne-par ticle Green's function in an external 
magnetic field p7| , [l5| . 

Consider first the case of an antiferromagnet. The per- 
turbation theory in the coupling constant for the parti- 
tion function (|ll]) can be developed in a standard way 
Jl5|-p7| . To this end we expand the square roots in a se- 
ries in the field it. The bare Green's function of this field 
has the form 
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G(°)(p,iuv) = -[c4 



+ pj + a(l - cosp 2 ) + / + h\ 



(15) 



where py = (p 2 , + P 2 ) 1 / 2 , w n = 2nn/u are the dimension- 
less Matsubara frequencies. In each order of perturbation 
theory one has to take into account all the possible con- 
nected diagrams with the fixed number of loops. 

We consider only the "renormalized classical" regime 



T> (max{/,a}) 1/2 c 



(16) 



since in the interval T < (max{/, a}) 1 / 2 c the staggered 
magnetization is well described by SSWT (sec a more 
detailed discussion of different temperature regimes in 
Ref. S). Since in this regime the quasimomentum cut- 
off parameter for quantum fluctuations is the boundary 
of the Brillouin zone A, while for thermal ones this is 
T/c, effects of quantum and thermal fluctuations can be 
separated (see e.g. Jl8|fl). Thus it is useful to perform 
the renormalization in two steps. At the first step the 
ground-state quantum renormalizations are performed, 
and the second step is the temperature renormalization. 
To this end we represent the renormalization constants 



Zi(g,u,A) = Z Qi (g,A)Zi(g r ,u r ) 



(17) 



where ZQi(g) contain ground-state quantum renormal- 
izatons, and Zi all the others, 



— Zqi9, u r — Z Qu U 



(18) 



are the quantum-renormalized coupling constant and di- 
mensionless inverse temperature. Note that the latter 
is renormalized due to the renormalization of the spin- 
wave velocity, c = Zq u cq. For further convenience we 
also introduce the quantum-renormalized anisotropy and 
interlayer coupling parameters 



f r — Zglf, a r — Z Q ^a 



Q3 



(19) 



which are just experimentally observed. Up to one-loop 
order we have 

Z Q = l-(N-l)^+0(g 2 ), 

Z Q1 = l-(N-2)^- + 0{g 2 ), Z Qu = l + 0(g 2 ) 

Z Q2 = l + ^ + 0(g 2 ), Z Q3 = 1 + ^ + 0{g 2 ) (20) 

Since the renormalization constants Zq^ are non- 
universal, the continual limit is insufficient to calculate 
them, and quantum-renormalized parameters can be de- 
termined only from the consideration of the original lat- 
tice partition function in the ground state. For 
the square-lattice antifcrromagnet this can be performed 
within the spin-wave theory, which is in fact a series ex- 
pansion in g (g ~ 1/5 for large 5). The results of the 



spin- wave theory are presented in Appendix A. We have 
to first-order in 1/5 



Z Q — V-^Qi — Z Q2 
Z Qu = 1 + 0.079/5. 



0.197/5, (21) 



Note that after performing the ground-state quantum 
renormalizations ( pf]| ) in the continual model (or, equiva- 
lently, renormalizations (|2l]) in the original lattice model) 
the theory becomes completely universal, i.e. thermody- 
namic properties do not depend on the cutoff parameter 
A. 

Now we pass to the consideration of the temperature 
renormalizations. The calculation of the renormalization 
constants Zj, is performed in the same way as in Ref. 
and leads to Z u = 1, 

Z =l+t r (N - l)ln(u r /z) 

+t 2 r (N - 1)(N - 3/2) ln 2 K^) + 0(4), 
Zi = l + t r (N - 2) \n(urn) + t 2 r (N - 2) ln 2 (w r /i) + Q(t 3 r ), 
Z 2 = l- 2t r ln(w r /i) + e>(i 2 ), Z 3 = l- U ln(u r fi) + 0{t 2 r ) 

where /it is an infrared cutoff parameter with the dimen- 
sionality of the inverse length and t r = g r / (2nu r ) . The 
only difference from the results of Ref. |L5) is that the 
ultraviolet cutoff parameter in ( |22| ) is u r (rather than 
A in the classical case), and two new renormalization 
constants, Z 2 and Z3, for the anisotropy and interlayer 
coupling parameters occur. 

The infinitesimal change of fj, generates the renormal- 
ization group transformation, and the derivatives of Z- 
factors with respect to /1 determine the renormalized- 
parameters flow functions (see, e.g., Ref. p7|). Since 
quantum-renormalization constants ( p(i| ) are invariant 
under RG transformation, it is sufficient to calculate the 
derivatives of Z. To the two-loop order we have 



f3(t r ) =^ = -(N - 2)t 2 r -{N- 2)t 3 r 

<;(t r )= f M^- = (N-l)t r + 0(t 3 r ) 

an 



0(t$) (23) 
(24) 



Other two flow functions will be needed only to the one- 
loop approximation: 



Ifitr) 

i^r ) 



/'- 



rfln Z 2 

d( 1 
d In Z 3 



dfi 



-2t r + 0{t 2 r ) 

-tr + 0(t 2 r ) 



(25) 
(26) 



Using (g3|)-(^6D we find the scaling laws for the Hamilto- 
nian parameters under RG transformation 



p = exp 



dt 
W) 



(22) 
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Z p = exp 



l/(JV-2) 



exp 



1 1 



N — 2 \t p t r 



[l+0(t„)]> 



where 
(27) 



0(t) 



eft 



AT- 1 
2(AT-2) 



(38) 



f p = f r exp 



atp = a r exp 



, x (N-l)/(N-2) 



If if) 



dt 



= fr\f 



-2/{N-2) 



is the sublattice-magnctization "critical exponent" in the 
temperature interval under consideration. This coincides 
wkh,,the e — ► limit of the critical exponent @2+e for 
d™2 + £ Q. The equation for t p /t r is given by @ 
which can be rewritten as 



[l + 0(t p) ], (29) k = l + trH k p N-2 ) 



UO(tp) 



(39) 



/?(*) 



dt 



x -l/(iV-2) 



where t p , f p , a p are the corresponding parameters of the 
effective Hamiltonian at the scale pp, and Z p is the scal- 
ing factor for the spin fields on this scale. 

Now we are ready to calculate the relative sublattice 
magnetization a r = ct/cto, where a = S/S (S = (Sq) is 
the staggered magnetization) and er = &{T = 0). The 
perturbation result in the zero magnetic field up to terms 
of the order of t 2 reads 

_ t r (N - 1) 2 
a r = 1 : m 



+ 



4 U 2A(/ r ,a r ) 
^(3-JVXJV-l), , 2 



32 

t 2 r (N - l)(B 2 - 2) 



ln z 



u 2 A(f r ,a r ) 



In 



u 2 A(f r ,a r ) 



(31) 



where 

" (32) 
(33) 

Note that the last term in (|3l]) corresponds to tem- 
perature renormalizations of the interlayer coupling and 
anisotropy parameters: 



A(f,a) = f + a +y/p + 2af 

B2=3 + f r /y/f? + 2a r f r 



ft = fr 



at = a r 



1 - t r In 



u 2 A(f r ,a r ) 



1 - -r In 



2 U 2 A(/rjar) 



(34) 
(35) 



We obtain from ( J28| ) the following scaling law for the 
sublattice magnetization 

a r (t r ,f r ,a r ,fi r ) = Zp 1/2 a r (t p ,f p ,a p ,np) (36) 

or, equivalcntly, 



(T r (tri fr: ®-ri AO 

Vr(tp,f P ,a p ,np) 



[l-P 2 (tr-tp)+0(t 2 p )] 



(37) 



Finally, the scale p is fixed by the condition that the 
■ar(g^ttients of logarithms in (31) on this scale are equal to 
unity, i.e. by a r (t p , f p ,a p , pp) = 1. Taking into account 
that u scales in a trivial way, u p — u r j p, we obtain the 
additional equation for p: 



2p 2 =u r A{f p ,a p ) 



(40) 



We have from @, (gg), © and @ the equation 
for the relative sublattice magnetization in the two-loop 
RG analysis 



W i/P2 =1 _ Z _L 
2 



+2(1 



(N - 2) In ■ 



u 2 A(f t ,a t ) 



jMl/ifr) (41) 

P2 



where / t and a t are the temperature-renormalized pa- 
rameters of the anisotropy and interlayer coupling. We 
derive from for these quantities 



ft/fr 

a t /a r 



-4/(AT-l) 



F 2/(iV-l) 



l + 0(t r /lT 1 r /' 32 ) 



(42) 
(43) 



The leading logarithmic term in the square brackets of 
( |4l| ) corresponds to SSWT, while other two terms de- 
scribe corrections to this theory. As it should be ex- 
pected, at low temperatures the RG results ([fl|), ( |42] ) 
and ( [f3[) reduce to corresponding perturbation expres- 
sions (|3l[), (HJ) and (|35|). We have to bear in mind that 

the quantity t r /or is a formal rather than real esti- 
mate for neglected terms since higher order terms also 
yield a contribution at not too low temperatures. Phys- 
ically, neglecting such terms is equivalent to neglecting 
3D (or Ising-like) fluctuations in the RG approach. 

As already pointed, the Neel temperature cannot be 
calculated directly in the RG approach since essentially 
non-spin-wave fluctuations contribute it and an accounce 
of diagramms with an arbitrary number of loops is re- 
quired. However, one can obtain a general expression 
for the Neel temperature in the following way. Consider 
first the temperature t* of the crossover to the true crit- 
ical region. This is determined by the condition t p ~ 1, 
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i.e. t* r-j a].^ 2 (the scale p is determined by (p 
stituting this into ( (4l| ) one obtains 

2 



Sub- 



t* = 2 



(N - 2) In 



u 2 r A{f*,a* 



21n(2/^) + C. 



AF 



(TV - 2) In -s— + 3 ln(2/f ncoi) + $af (0) 
u£a r 

(51) 



where /* = f r (t* 



or 



a r (t*) 



l/(N-2) 



(44) 
Caf is 



(hi) large- iV case. To retain the structure of the equation 
(E) at finite N, it is convenient to expand (^0|) in 1/(N— 
2jrather than in l/N. To first order we obtain 



some constant of order of unity. With further flow of 
RG parameters, 3D Heisenberg (or 2D Ising) fluctuations 
change only the constant Caf which is replaced by the 
universal function $AF(a r // r ) ~ 1. Thus for the Neel 
temperature we have 



A(/ r ,a r ) 
2 

A(/ r ,a r ) 



Bo — 2 , t r 
1 + — In — 

N -2 t„ 



(B 2 -2)/(JV-2) 



(52) 



where 



(N - 2) In 



u 2 A(f c ,a c ) 



21n(2/t Neel ) 



Theip. we have 

*AF(a r //r) 



f - f i 2/(Ar ~ 2) a 



a t 1/{N ~ 2) 



(45) 



(46) 



[recall that t = T/(2?rp s ), u r = c/T]. The function $ 
is determined by non-spin-wave fluctuations and cannot 
be calculated within the RG approach. Note that the 
equation (^H|) is analogous to the result for the correlation 
l ength in the renormalized classical regime in the 2D case 
|p.7|,p.8| , which may be rewritten in the same manner: 

tr = 2 / [(N - 2) H£ 2 /u 2 r ) + 2 ln(2/t r ) + In Q] (47) 

where is an universal numerical preexponential factor. 
Comparing (|4^) and ([17]) we see that the only difference 
is that in the quasi-2D anisotropic case the logarithms are 
cut at A 1 / 2 (rather than at l/£ in the isotropic 2D case). 
In the quasi-2D isotropic case the quantity $ can be be 
calculated within the l/N expansion (see the next Sec- 
tion) , and a more general case requires numerical analysis 
(e.g., the quantum Monte-Carlo method). 

Three limiting cases may be considered: 
(i) anisotropic 2D case a — 0. The equation for magneti- 
zation (p4j) takes the form 



2 



(JV-2)ln-l- + -iln(l/ff P ) 

Urfr P2 



(48) 



+2(l-a 1 J^) + 0{t r /a 1 r l^) 

The equation for the Neel temperature (|45|) reads 

1 



^Ncol — 2 



(N - 2) In - jy + 4 ln(2/t Neel ) + $af(oo) 

(49) 



U 2 /r 



(ii) isotropic quasi-2D case / = 0. We obtain 



gi/ft = 1 _ i 
2 



(N - 2) In ■ 



-ln(lM.) 

P2 



(50) 



+2(i-ay' 32 )+o(t I ./ay' 32 ) 



and 

*Neel 



1 



(iV - 2) In ■ 



+2(1 



U 2 A(/r,«r) 

a 2 r ) + 0{t r /a 2 r ,l/N) 



B 2 ln(l/a 2 ) (53) 



(JV - 2) In 



u 2 A(/ r ,a r ) 



B 2 ln(2/t 



Nccl 



+ $AF((* r //r) 
(54) 



The ferromagnetic case can be considered in a similar 
way. In this case the quantum ground-state renormaliza- 
tions are absent and the "classical" regime (an analogue 
of the "renormalized classical" one in the antiferromag- 
netic case) is determined by 



T> JSm&x{f,a} 



(55) 



To develop perturbation theory for the partition func- 
tion (^) we pass from the real fields ir x , n y to the cyclic 
components 



7T r ± ITT,, 



(56) 



and expand again square roots in n + , tt . The bare 
Green's function of the fields ir + , ir~ takes the form 



G (0) (p, iuj n ) = - iijj n +p\ + a(l - cospz) + / + h 



(57) 



Due to absence of quantum renormalizations we have 
Zqi = 1 and the indices r may be dropped. The factors Zi 
has the same form ( p2] ) as in the antiferromagnetic case 
with N = 3 and the replacement ln(u r /i) — > ln(u/i)/2. 

The relative magnetization a = S/S (S = (S z )) to the 
two-loop approximation reads 



1 



-1 2 

2 n uA(f,a) 



-(52 -2) In 



uA(f,a) 



(58) 
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(t = g/u). The scaling equation (f37|) is valid in the fer- 
romagnetic case too. The equation for p in the this case 
takes the form 



2p 2 = uA{f p ,a p ) 



(59) 



Thus we obtain the equation for the magnetization in the 
form 



(7 = 1 — 



t 



In- 



uA(f t ,a t ) 



21n(l/a) 



(60) 



+2(1 - 5) + 0{t/a) 



The results for the temperature renormalization of the 
anisotropy and interlayer coupling parameters have the 
same form (fl2"|), (fl3|). The Curie temperature is deter- 
mined in the same way as the Neel temperature in the 
antiferromagnetic case. The result reads 



fa 



= 2 



In 



uA(f c , a c ) 



+ 21n(2/t Cu ric) + *F(a//) 



(61) 



The calculation of the magnetization and ordering tem- 
perature of a classical magnet is performed in Appendix 
B. 

Thus the RG approach is sufficient to calculate the 
magnetization in the spin-wave and the 2D-like regions 
and to calculate the Neel (Curie) temperature up to some 
universal constant. The crossover temperature region of 
the quantum antiferromagnet can be considered within 
the 1/N expansion. Besides that, in the case of quantum 
quasi-2D antiferromagnet, this expnasion enables one to 
describe the true critical region and to evaluate the quan- 
tity $af(0). 



IV. COMPARISON WITH THE 1/N EXPANSION 
IN THE QUANTUM O(N) MODEL AND THE 
CROSSOVER TO THE CRITICAL REGIME 

The 1/N expansion gives a possibility to develop an- 
other perturbation theory for the partition function (11). 
Unlike the renormalization group approach, this method 
works satisfactorily at arbitrary temperatures. As well 
as in Section 3, we consider only the ordered phase. 

Consider first the case of the antiferromagnet. We be- 
gin with the generalization of the results of Ref. || to 
the case where anisotropy is present. To develop pertur- 
bation theory in 1 /N we integrate out a- fields from (11). 
Thus we have 



Zai 



D\exp{NS e ff[\,h}) 



Seff[\h] = ilndetG + ^(l-a 2 )Sp(zA) 

+ ^S P [(iA5 - h/p*) Go - h/p°) 



(62) 



(63) 



where 

G™" 1 = \d 2 T + V 2 + aA z + /(l - S^)}- 1 (64) 
A z a lz (r,r) = ir^+i (r,r) - a iz (r,r) 

and 5 = (<7jv(r, r)) is the relative staggered magnetiza- 
tion. 

Since N enters ([32]) only as a prefactor in the exponent, 
expanding near the saddle point generates a series in 1 /N. 
To zeroth order in 1/N the excitation spectrum, which is 
given by the poles of the unperturbed longitudinal and 
transverse Green's functions, contains a gap f 1 / 2 for all 
the components a m except for m = N: 

-l 



G?(q,w n ) 



4 + q 2 + a(l-COSq z ) + f 



q 



a(l — cos q z 



(65) 
(66) 



The absence of the gap for the N-th (longitudinal) mode 
in the ordered phase is an exact property of the model 
under consideration in any order in 1/N. 

The sublattice magnetization at T < Tn cc i is deter- 
mined by the constraint equation (cr 2 ) = 1. To first order 
in 1 /N the constraint takes the form 

! -2 N-l ^ f d 2 k|| j dk z 



(27T) 2 

"* — 7T 

+g [F(T,a)-R(T,a)], 



(67) 



where 



7T 

R(T^ = T^J^J^[GU^ m )] 2 (68) 



[£t(k,w m )-E,(O,0)] 



d 2 k|. 



2tt n mj n(k, Wm ) 



(69) 



The longitudinal and transverse mode self-energies are 
given by 



£ 4) i(k,w m ) 



2T^ 
~N~ ^ 



d\ \ f dq z G° t l (k-q, 

(27T) 



2tt 



n(q,^„) 



where 

n(q,w„) = II(q, w n ) 



-=2 



2a 



■G?(q,w n ), 



(70) 

(71) 
(72) 



N - 1 



N 



G°(p, W/ )G°(q + p,^, 



+-jyG?(p,wj)Gj , (q + p,w n + wi) 
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Since the polarization operator Tl(q,uj n ) enters only the 
first-order corrections in ( |67| ) , the contribution from lon- 
gitudinal Green's functions to II influences thermody- 
namic quantities in order of 1/N 2 and can be formally 
neglected. However, one should bear in mind that for 
finite N this corrections may be large. Physically, the 
neglection of the longitudinal part of IT corresponds to 
neglecting the contribution of Ising-like (spin-flip) exci- 
tations. However, these excitations are negligible only at 
wavevectors q J 1 / 2 (see, e.g., Ref. Pq|). As follows 
from (|72|), such quasimomenta yield a dominant contri- 
bution at T < 2a 2 [(a + f)/f} 1/2 /g. The opposite case 
T 3> 2a 2 [(a + /) / /J 1 / 2 / g corresponds to the Ising critical 
region which cannot be treated within the 1/N expan- 
sion. 

As well as in the previous Section, we consider only the 
case T ^> (max{/, a}) 1 ^ 2 c. The procedure of integration 
and frequency summation in ( |67| ) is analogous to that of 
Refs. [p)J|. We obtain 



1 - 
-2 



T 



(„- 2)1 „^! + S2 ,„MTV^A) + lT 
cr A Xtf 



Airp 
ln(2T 2 /c 2 A) 



ln(2T 2 /c 2 A) + x- 

1 ln(2T 2 /c 2 A) + x w 

! + T7 m h&? 

N Xtf 



(73) 



where A = A(f r ,a r ), (see Eq. (|32j)), B2 is determined 
by®, 



^Ps _ 2 



(N — 2)T 



(74) 



and we have introduced the quantum-renormalized pa- 
rameters 



f r = /(l - 2Q t 



a(l - Qf 



(75) 



(l+4Q A )p; 



N=oo 



gp s {l-Qiv)/N (76) 



with =°° = Nc(l/g - A/2tt 2 ) 
ized spin stiffness in zeroth order 



is the renormal- 
in 1/N, Qa = 
(8/37r 2 iV)ln(A^Ac/16p s ) (in this Section we use the rel- 
ativistic cutoff uj 2 + k 2 < A 2 of frequency summations 
and quasimomentum integrations; for this regularization 
scheme the bare spin wave velocity Co is replaced by the 
quantum-renormalized one, c). Since another regulariza- 
tion scheme is used, the expressions (|75|), ( f76|) are differ- 
ent from the corresponding results of Sect. 3. As well as 
in Sect. 3, the quantum-renormalized parameters are not 
universal and therefore should be determined from the 
spin- wave theory (see Appendix A) rather than from the 
continual model. The functions Ii,2(x) are some func- 
tions with the asymptotics 1/x at large x, so that at 
3> 1 their contributions are small. For the isotropic 
quasi-2D case these functions were calculated in Ref. [|| . 
Consider first the case of not too high temperatures 



T{N - 2)/Anp s < a 2 . 



(77) 



where Ii^(xw) are small enough. Using the identity 
ln(T 2 /c 2 A) + x w = Att Ps /(N - 2)T which holds to ze- 
roth order in 1 /N, we transform the logarithmic term in 
the right-hand side of (|7^) into a power to obtain 



[l-h(xff)]a 



(N-2)ln^ 2 



1 - 



T 



2T 2 



Airp s 
B 2 ln(l/a 2 ) 



2 + 2a 2 - h{x w ) 



(78) 



Note that in the quasi-2D case (/ = 0) the equation ( JTS 
slightly differs from the equation (54) of Ref. |3j by the 
replacement 



1-S 2 



ln(2T 2 /a r )/[ln(2T 2 /a r ) + x„] 

which holds to zeroth order in 1/N. 

At temperatures T(N — 2)/Anp s <C a 2 the contri- 
bution of I\ , I2 can be neglected and the fluctuations 
in this temperature region have a 2D-like Heisenberg 
nature. In particular, in the low-temperature region 
T(N - 2)\n(2T 2 /A)/4:TTp s < cr 2 the result of SSWT 

ran 



T(JV-l), 2T 2 

= 1 " In 

87rp s c 2 A 



(79) 



is reproduced. One can see also that at I\ = I2 = the 
resul t (|78| ) coincides with the large- N limit of the RG re- 
sult (|4j) (see Eq. (||)). However, at finite N the renor- 
malization group provides a more correct description of 
the sublattice magnetization at T(N — 2)/Anp s <C a 2 . 

In the temperature region a 2 ~ T(N — 2)/Aitp s , which 
corresponds to the crossover to the true critical behav- 
ior, the situation changes. In this case the large-iV re- 
sult demonstrates a more sharp decrease of a than the 
RG approach. In the quasi-2D case || the result (|7^) is 
smoothly joined with the 3D temperature dependence 
(see below). Thus in the quasi-2D case the result of 
the 1/N expansion should be considered as an interpo- 
lation between the low-temperature and critical regions. 
One could expect that this holds also in the presence of 
the anisotropy were critical behavior cannot be described 
within the 1/N expansion. 

In the region < 1, i.e. ct 2 < T(N-2)/4np s the true 
critical behavior takes place. In the isotropic quasi-2D 
case (/ = 0) the result for the staggered magnetization 
reads (see Ref. pf) 



ATrp s 



(N - 2)T No ci 



Pa /Pa 



1 



1-Aq 



T 



2ih 



(80) 



where (3% = (l — 8/w 2 N^ /2 is the true 3D critical ex- 
ponent for the order parameter (for N — 3 we have 
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0a ~ 0.36), A = 
pansion for Tn co i 
with 



2.8906/A. The result of the 1/N ex 
H satisfies the general formula 



$af(0) = -0.0660 



(81) 



In the anisotropic 2D case (a — 0) the temperature de- 
pendence of magnetization in the critical region is deter- 
mined by the Ising-like excitations (domain walls) which 
cannot be considered within the 1/N expansion. The 
universality hypothesis predicts the same temperature 
behavior of a as in Ising systems 



A(l - r/T Nee l), 



(82) 



where A is some constant. As demonstrate numerical 
calculations (see the next Section), the temperature de- 
pendence of a determined from ( |T8|) is smoothly joined 
with d82| ) , A and Tn cc i being considered as fitting param- 
eters. 

In the presence of both anisotropy and interlayer cou- 
pling, the situation in the critical region x-& < 1 is more 
complicated. Consider first the case / < a. Then at 
1 > Xa > {f /(a + Z)] 1 / 2 the (sublattice) magnetization 
has the 3D Heisenberg behavior ( pC| ) with some coeffi- 
cient 1/[1 - A (f/a)]; at x w < [f/(a+ /)] 1/2 the behav- 
ior a(T) changes to the 3D Ising one. At / ~ a the 3D 
Heisenberg region disappears, and in the whole critical 
region the 3D Ising behavior takes place. With further 
increase of / (at / > a), the 2D Ising critical region oc- 
curs for 1 > x-& > Xo{f /a), while at x-^ < xo(f /a) the 3D 
Ising behavior still takes place. However, the dependence 
Xo(f /a) cannot be calculated within the approaches un- 
der consideration. 

Now we turn to the case of a ferromagnet. As already 
mentioned, in this case the dynamical part of the action 
cannot be generalized to arbitrary N. Thus the expres- 
sions (|67|)-(|72l) of the first order in 1/N (where we put 
N = 3) should be considered as physically reasonable 
rather than strict results. The dynamical part of the 
action in (^) results in cutting the quasimomentum inte- 
grals at q ~ {T/JSf/ 2 . The indices r can be dropped, 
since the quantum renormalizations are absent. Thus we 
have instead of (78) for T ^> JSmax{f, a} 



obtained in the practically interesting case N = 3. In the 
spin-wave and 2D-like regions, i.e. at 



o> > T/4wp s , T > A 



(84) 



we have the RG result for the relative (sublattice) mag- 
netization 



T 
4irp s 



2T(T) 

In , , „ V \ + 2m(l/o>) + 2(1 - o>) 



A(f t ,a t ) 



(85) 



where the function A(/, a) is determined by (p2[), the 
temperature-renormalized values of interlayer coupling 
and anisotropy parameters are 



ft/fr = (at/otr) 



2 -w 2 



(86) 



and the quantities T(T),a r , f r ,ot r , p s are given in the Ta- 
ble 1 (see also Appendix A). 



T(T) 0Y p s f r a r 
quantum AFM T 2 /c 2 S_/S Q jSS fsl/S 2 aS /S 
quantum FM T/JS S/S p Q s f a 

classical FM, AFM 32 S/S p° s Z L1 fZ^ aZ£^ 



Corresponding equation for the magnetic ordering tem- 
perature Tm has the form 



T 



M 



hi 



2T(T 1 



M 



A(/ c ,Oc) 



2 In 



^Ps 



Hf/a) 



(87) 



where $(x) is some function of order of unity (in the 
quantum case it is universal, i.e. does not depend on the 
upper cutoff parameter) , f c and a c are the temperature- 
renormalized interlayer coupling and anisotropy param- 
eters at T = Tm that are determined by 



fc/fr = (a c /a r ) 2 = (Tm/^Ps) 2 



(88) 



a=l- 



T 
4irp° s 



In- 



2T 
JSA 



B 2 HI/ a 2 ) - 2 



2a 2 



0(47rp° s 
(83) 

where 0(Aitp1/lj 2 ) terms cannot be calculated within 
such a consideration. For the Curie temperature we re- 
produce the RG result (|6l|). 



V. DISCUSSION AND COMPARISON WITH 
EXPERIMENTAL DATA 

The above consideration provides a description of the 
long-range order of quantum and classical magnets in dif- 
ferent temperature regions. Let us summarize the results 



Since Tm/^Ps ~ 1/ m (l/A) <C 1 the temperature renor- 
malizations are important when treating experimental 
;_da|ta. In particular, the parameters, which are measured 
' 17 ax, different temperatures, may differ considerably. 
In the case a = we have 



4irp s 



T 



M 



Anp s 



hi 



ln^2^+41n(l/ci r 
Jr 

r(T M ) 



fr 



41n 



47T / O s 

T M 



2(1 -o> 



*(0) 



In the case / = we obtain 
T 



a r = 1 



4irp s 



In Si + 31n(l/cr r ) + 2(1 - a r ) 



T M = 4:Trp s 



ln^M + 31n^ + (f >M 

a r Tm 



(89) 
(90) 

(91) 
(92) 
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The results of solving the SSWT equations |1||2[ differ 
from ©-(H) by the replacement 4(3) -> 2(1) for the co- 
efficient at the second term in the square brackets (which 
yields the double-logarithmic correction to the standard 
SWT) in the anisotropic 2D (isotropic quasi-2D) case re- 
spectively. Thus the role of the corrections to SSWT is 
more important in the isotropic quasi-2D case than in 
the 2D anisotropic one. 

As discussed in the Introduction, the results (|B|), j87| ) 
are expected to hold in the first order of the 1 /M expan- 
sion in the CP* 1 ' 1 model (at M = 2). 

In the temperature interval outside the critical region 



a 2 r > T/4np s , T> A 



(93) 



the result of the 1/iV expansion in the O(N) model to 
the first order in 1/N reads 



[l-h(x w )]a r = 1 
+2B 2 ln(l/ov) + 2(1 



T 
4irp s 
a 2 ) 



In 



2r(T) 



A(/ r ,a r ) 
h(x<r)] (94) 

where x-& — Airp s a 2 /T, B 2 and A are determined by ( |33[ ) , 
([32|), Ii, 2 (x) are some functions with the asymptotics 1/x 
at large x, other quantities are given in the Table 1. In 
particular cases a = and / = the coefficient at the 
second term in the square brackets in (|94| ) is too times 
larger than for the RG results (|S9"|), (|9l[). In the spin- 
wave and 2D-like temperature regions this is an artifact 
of the first-order 1/N expansion. At the same time, the 
1 /N expansion provides a more correct description of the 
crossover temperature region. Due to the difference in 



the crossover conditions (93) and (84), the equations for 
Tm have the same form (30), ( p2] ) in both approaches. 

Now we discuss the experimental situation. First we 
consider the temperature dependence of the sublattice 
magnetization in La2Cu04 (Ref. [Q) which is shown in 
Fig.l. This Figure presents also the results of spin- wave 
approximations (SWT, SSWT and the Tyabhkov theory 
Ipjfl ; see a more detailed discussion in Ref. ||), RG ap- 
proach and the result of 1/iV expansion (|9~4|). The value 
7 ~ 1850i^ was calculated by using the experimental 
data (3^] while a r = 1-10 -3 was obtained from the best fit 
of experimental dependence a r (T) to the spin- wave the- 
ory at low temperatures. The result of the 1/N expansion 
to first order in 1/N is 7n cc i = 345 K which is consid- 
erably lower than for all the spin-wave approximations 
and is in a good agreement with the experimental value, 
^Noci = 325K. The RG approach describes correctly the 
dependence o>(T) in the spin- wave region (T < 3001^) 
and 2D-like region (which is very narrow since a is very 
small) while at higher temperatures this approach over- 
estimates a. At the same time, the 1 /N expansion curve 
is most close to the experimental data and demonstrates 
a correct critical behavior. The results of the numerical 
solution of equation ( |94|) in the temperature region ( p^ ) 
and the dependence (BOh in the critical region turn out 
to be smoothly joined at the point T — 330K (marked 
by a cross). 



In the crossover region (320K< T < 340K) the theoret- 
ical 0(3) curve lies slightly higher than the experimental 
one. One may speculate that this is due to the influence 
of anisotropy. Fixing A in (94) and determining B 2 from 
the best fit at intermediate temperatures (see Fig.l) one 
finds the values a r = 1 ■ 10~ 4 , f r = 5 • 10~ 4 . This value 
of a is more close to the experimental data of Ref. |52) . 
Thus our approach gives a possibilty to estimate the rel- 
ative role of interlayer coupling and magnetic anisotropy 
in layered compounds. 

In the layered perovskites K2MF4, Rb2NiF4 and 
K2MnF4 the magnetic anisotropy is known to be more 
important than the interlayer coupling. K2MF4 has spin 
5 = 1, and neutron scattering data yield \J\ = 102K and 
^Nooi = 97. IK (see Ref. U). Fig. 2 shows the experimen- 
tal dependence <r(T) |^6| and the results of the spin-wave 
approaches, RG approach and the numerical solution of 
Eq.@. The value f r = 0.0088 was obtained from the 
best fit of the result of SSWT to experimental data at low 
temperatures (this value coincides well with the experi- 
mental one f r = 0.0084, Ref. M). In the spin-wave and 
2D-like temperature intervals (§3) (T < 80K) the curves 
corresponding to the 1/N expansion and RG approach 
lie somewhat higher than the experimental points since 
T 2 / f r c 2 in this region is not large, and the renormalized- 
classical description is not too good (a more accurate cal- 
culation can be performed by carrying out exact summa- 
tion over the Matsubara frequencies). Bearing in mind 
this correction, the RG approach gives a more correct 
qualitative tendency than the 1/iV-expansion in the 2D- 
like region. At the same time, the 1/N expansion curve is 
in a good numerical agreement with experimental data. 
The joining procedure with the Ising critical behavior 
(p2|) may be performed in a rather wide temperature 
region 0.85T Nccl < T < 0.9T Nccl and gives A = 0.01, 
Tkcci = 91. 4K. The width of the critical "Ising" region 
makes up about IK. Note that an account of the terms of 
order of 1/x^ in ( |78| ) , which can be performed by analogy 
with the calculations of Ref. ||, gives Tn oc i = 92.7K. 

In the crossover region (80K< T < 90K) the the- 
oretical 0(3) curve for K2MF4 lies, in contrast with 
the case of La2Cu04, slightly lower than the experimen- 
tal one. This fact may be attributed to the influence 
of interlayer coupling. The fitting in the crossover re- 
gion yields the values a r — 0.0017, f r — 0.0069 which 
correspond to TnccI = 97K and the bare parameters 
a\J\ = 0.1K, £|J| = 0.76K. Direct experimental data 
for a are absent, but our estimation seems to be reason- 
able. 

Rb2NiF4 has a larger magnetic anisotropy. According 
to Ref. §J, one has |J| = 82K, \ J\f r = 3.45K, = 
94. 5K. From the best fit of SSWT to the dependence 
a r (T) at low temperatures one obtains f r — 0.046 which 
is also in a good agreement with the above experimental 
value. Then one obtains from ( p4] ) Tn cc i = 95. 5K. 

K 2 MnF4 has spin S = 5/2 and therefore represents 
a situation which is intermediate between the quantum 
and classical cases. Fig. 3 shows a comparison of the re- 
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suits of different approaches with experimental data for 
this compound. The parameters used are |J| = 8.4K, 
\J\f r = 0.13K (see Ref. 0). One can see that the l/N ex- 
pansion yields good results, and the experimental points 
lie between the quantum and classical RG curves, the 
quantum approximation being essentially more satisfac- 
tory. This confirms once more that it is difficult to realize 
the classical limit (see Appendix B). Note that SSWT, 
which correctly takes into account lattice effects, provides 
in this case better results in comparison with the RG ap- 
proach. Thus an accurate treatment of such situations 
within continual models requires numerical calculations 
of quasimomentum integrals and sums over Matsubara 
frequencies in (67). 

Fig. 4 shows a comparison of the results of SSWT and 
the RG approach for the magnetization of a classical mag- 
net with the Monte-Carlo calculations jl3| . One can see 
that, except for a very narrow critical region, the RG 
curve is rather accurate although topological excitations 
are neglected. Note that the region of applicability of the 
RG approach in the classical case is more broad than in 
the quantum case, so that we need not to use the large- 
N approach for describing the crossover to the critical 
region. 

Thus the RG approach (or, equivalently, the l/M ex- 
pansion in the SU(M) model) and the l/N expansion in 
the O(N) model turn out to give good results in different 
temperature regions. Whereas the first method describes 
well the 2D-like regime, the l/N expansion describes suc- 
cessfully the crossover to the critical region and, in the 
absence of anisotropy, the critical behavior. Both meth- 
ods give the same results for the ordering point up to 
the terms of order of lnln(l/A). Besides that, the l/N 
expansion permits to calculate nonsingular terms in the 
quasi-2D case. 

To conclude, our results give a possibility to describe 
magnetic properties of real layered magnets with a rather 
high accuracy. The approaches applied may be useful 
for treating magnetic and structural phase transitions in 
systems with more complicated order parameters |33|] . 



p s =jSS , c 



V&yS 



(96) 



with 7 being the quantum-renormalized intralayer ex- 
change parameter determined by 



l/\J\ 



1 



-y 

k 



1 



1 



1 + ^ (97) 



For the quantum-renormalized coupling constant we have 
g r = c/p s . The quantum-renormalized interlayer cou- 
pling and anisotropy parameters can be determined from 
the first-order l/S'-corrections to the excitation spec- 
trum. We have at T = and for small in-plane wavevec- 
tor components E2[ 



El 



8( 7 5) s 



2V 



1 — cos fe. 



_6_ 



where 



7 



%(S /S)\J\ 



is the renormalized interlayer coupling and 
6=(So/S) 2 {(2S~l)( + 47 1 S\J\/ 1 } | J | 



(98) 



(99) 



(100) 



is the renormalized anisotropy. Note that in the case 
£, T] <C 1, which is considered only, single- and two-site 
anisotropies lead to the same effects. Comparing the 
spectrum ( B8|) w ith the bare spin-wave spectrum deter- 
mined froniflq) 



El=c 2 



k 2 + a(l - cos k z ) + f 



(101) 



we obtain the relation between the bare and the 
quantum-renormalized parameters : 



fr = -s = (So/S) 2 



(2S-1)C|J| , A V J 2 



lS 



T 



a r 



2i 



aSo/S 



(102) 
(103) 



APPENDIX A. SPIN- WAVE RESULTS FOR THE 
GROUND-STATE RENORMALIZATIONS IN A 
QUANTUM ANTIFERROMAGNET. 

The ground-state thermodynamic quantities of the 
quantum antiferromagnet can be calculated within the 
spin-wave theory. The result for the ground-state stag- 
gered magnetization reads (see, e.g., ||) 



s ° - s ~ 2 y 

k 



1 



S- 0.1966 (95) 



where </>k = ^(cosk x + cos fey). The ground-state spin 
stiffness and spin- wave velocity to first order in 1/S are 
given by g 



However, it should be noted that since the spectrum ( |9q ) 
contains only renormalized parameters rather than the 
bare ones, (a/2)|J| and [2£(1 - 1/25) + Arj\ \ J\, only 7' 
and S can be determined experimentally. 



APPENDIX B. RENORMALIZATION GROUP 
ANALYSIS IN THE LATTICE O(N) MODEL AND 
THE LIMIT OF CLASSICAL SPINS 

The treatment of the partition function for the classical 
anisotropic quasi-2D anisotropic magnets (|^) is similar 
to the isotropic 2D case In this case the relative 

temperature t = Tj (2irp° s jplays the role of a coupling 
constant, and we have instead of the first line of (|l^) the 
scaling relation 



11 



t — Z\Ir 



(104) 



where tn is the renormalized temperature. The bare 
Green's function of the field it = n (nz)z has the form 



^(q) = ^[2(2- 



Note that in the case a = a similar expression was ob- 
tained in Ref. |p2| . However, the result of this work con- 
tains wrong coef ficient at the second term in the square 
brackets of (113) since not all two- loop corrections were 
taken into account. 



cosg^a — cos q y a) 
i-i 



+a(l - cos q z a) + f + h]~ 



(105) 
(106) 



where a is the lattice constant. The renormalization con- 
stants can be calculated from the two-point vertex func- 
tion. It is useful to represented this constants as 



Comparing the result for the magnetization (112) with 
these of Sect.2, (|l|) and ©, we can write down the 
criteria of applicability of the classical limit (see Sect.2) 
with the correct numerical factors 



T 2 > 32c 2 , 
T > 32 JS, 



AFM, 
FM. 



(114) 



Zi(t,a) = Z Li (t)Zi(t L ,a) 



(107) 



where = tZ Ll , Zu contain non-logarithmic terms 
which are not changed under RG transformations, and 
Zi contain all the other terms. We have 



Zl 



Zli 
1 



Z L3 = l- 7rf/2 + 0(i 2 ) 



(108) 



The results for Zi read 

Z = l + t L (iV-l)ln(64a/i) 

+t 2 L (N - 1)(JV - 3/2) ln 2 (64a/j) + 0(t\), 
1 + t L {N - 2) ln(64a/i) + t\ (N - 2) In 2 (64a//) 



Zi 



It is difficult to satisfy these criteria in the ordered phase 
T < T M ~ l/ln(l/A) at not too small A and 1/5 due 
to the large value of the numerical factor in (114). 
Figure captions 

Fig.l. The theoretical temperature dependences of the 
relative sublattice magnetization <x r from different spin- 
wave approximations, RG approach (Eq. (|9l|)) and 1/N 
expansion in the O(N) model (Eqs. ( fT8[ ) and (|80|)), and 
the experimental points for La2Cu04 (Ref. ]3(J). The RG 
curve is shown up to the temperature where the deriva- 
tive da r /dT diverges. The curve denoted by 1/N' is 
the best fit in the crossover temperature region to the 
experimental data with the anisotropy being the fitting 
parameter (see discussion in the text), 
ij Fig. 2. The relative staggered magnetization a r (T) 



Z 2 = l- 2t L ln(64a/i) + 0(t 2 L ), Z 3 = l-t L ln(64o M ) + 0[t L J, ^iffif 4 (points) as compared to the standard spin 



The expression for the magnetization to two-loop approx- 
imation is 



a = 1- 



t L (N-l) 



In - 



G4 



*i(3- 



N){N 



-1) 



t 2 (N 



32 
1) 



In 2 . 



G4 



A(f L ,a L ) 
h 



y/fl + 2a L f L 



In- 



64 



A(/l,«l) 



where we have defined the renormalized quantities in the 
lattice case, 



/l = fZ L 2 , a L =aZ 



L3 



(111) 



The equation for the magnetization can be derived in the 
same way as in Sect.2 to obtain 



<fV/3 2 = 1 _ l ± 

2 



(N - 2) In 



+2(1 



G4 



-ln(l/a) 

P2 



A(/ t ,a t 



(112) 



wave theory (long-dashed line), SSWT (dot-dashed line), 
RG approach and result of the solution of ([78]) in the 
intermediate-temperature region ( |93|) (solid line). The 
short-dashed line shows the extrapolation of the 1 /TV- 
expansion result to the Ising-like critical region according 
to Eq. (|32"|) . The boundary of the 2D-like and crossover 
regions is marked by an arrow. 

Fig. 3. The dependence a r {T) for K2MJ1F4 (points) as 
compared to the results of SSWT (dashed line), RG anal- 
ysis in the quantum (dot-dot-dashed line) and classical 
(dot-dashed line) limits and solution of ([78|) (solid line). 

Fig. 4. The renormalization group (solid line) and 
SSWT (dashed line) results for the relative magnetization 
a of a classical anisotropic 2D magnet (£ — 0, 77 = 0.001) 
as compared to the result of the Monte-Carlo calculation 
U|. The RG and SSWT curves are shown up to the 
temperature where da/dT = 00. 
Table caption 

Table 1. Parameters of the equations for the (sublat- 
tice) magnetization (|85|), ( |94| ) for different cases, Zli = 
Z L 2 =Z L3 = 1- T/8np° s . 



Then the ordering temperature satisfies the equation 

64 



tM = 2 



(N - 2) In 



A(/ c ,Oc) 



+ 21n(l/< M ) + $d(a//) 
(113) 
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